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Abstract 

We present the analytic next-to-next-to-leading perturbative QCD corrections in the leading 
twist approximation for the moments N = 2, 4, 6, 8 of the flavour singlet deep inelastic structure 
functions F 2 and F L . We calculate the three- loop anomalous dimensions of the corresponding 
singlet operators and the three- loop coefficient functions of the structure functions Fl and F 2 . 
In addition, we obtained the 10 th moment for the non-singlet structure functions in the same 
order of perturbative QCD. We perform an analysis of the obtained results. 

1 Introduction 

The calculation of the next-to-next-to-leading (NNL) QCD approximation for the structure func- 
tions i*2 and Fl of deep inelastic electron-nucleon scattering is important for the understanding 
of perturbative QCD and for an accurate comparison of perturbative QCD with experiment. To 
obtain the NNL approximation for these structure functions in the operator product expansion 
(OPE) formalism one needs the 3doop anomalous dimensions of the operators, the 2doop Wilson 
coefficient functions for Fi and the 3doop coefficient functions for Fl. At present, these struc- 
ture functions are known in the next-to-leading approximation only, since the 3-loop anomalous 
dimensions and the 3-loop coefficient functions for Fl were not calculated yet. 

The 1-loop anomalous dimensions were calculated in Ref. |l]]. The complete 1-loop coefficient 
functions were obtained in Ref. || (see also the references therein). Anomalous dimensions in 
2-loop order were obtained in Refs. ||-||] and the 2-loop coefficient functions were calculated in 
Refs. [ffl-[|ll||. In a previous paper [Q| we presented the NNL corrections of the non-singlet type 
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in the leading twist approximation for the moments N = 2, 4, 6, 8 of the deep inelastic structure 
functions F<i and Fl- 

In the present paper we calculate the NNL QCD corrections to the singlet moments N = 2, 4, 6, 8 
of both structure functions F2 and Fl. To this end, we calculate the corresponding 3-loop anomalous 
dimensions and the 3-loop coefficient functions for the structure function Fl- In addition, we 
present the 3-loop coefficient functions for the structure function F2 for N = 2,4, 6, 8. We also 
obtain the N=10 non-singlet moments of F2 and Fl- The calculations are done for the leading 
twist approximation for zero quark masses. 

In the fifth part of this paper we analize the effects that the calculated 3-loop corrections have 
on the structure functions. 



2 Formalism 

We need to calculate the hadronic part of the amplitude for unpolarized deep inelastic electron- 
nucleon scattering which is given by the hadronic tensor 

W^(p,q) = ^|d 4 ze^(p,nucl|J M (z)J,(0)| nucl.p) 

= e^F L (x,Q 2 ) + d^—F 2 (x,Q 2 ) (1) 



(Ax 2 2x \ 

where is the electromagnetic quark current, x = Q 2 /(2p ■ q) is the Bjorken scaling variable 
(0 < x < 1), Q 2 = — q 2 is the transferred momentum and | nucl,p) is the nucleon state with 
momentum p. Spin averaging is assumed. The longitudinal structure function Fl is related to the 
structure function F% by Fl = F2 — 2xF%. 

As one approaches the Bjorken limit, Q 2 — > 00, x fixed, one can show that the integration 
region in Eq. (||) near the light cone z 2 ~ progressively dominates fl3]| , due to increasingly rapid 
phase fluctuations of the term e tq ' z outside the light cone region (and presuming that the integrand 
(p, nucl| Ju(z) e/jy (0) I nucl,p) varies smoothly outside the light cone). Since we have to deal with this 
non-local limit z 2 ~ 0, a formal operator product expansion in terms of local operators can only 
be applied together with the dispersion relation technique [14|. These techniques together provide 
a systematic way to studyf] the leading and non- leading contributions to the hadronic tensor. 

The tensor Wu V is, by application of the optical theorem, related to a scattering amplitude T^ u 
which is a more convenient quantity for practical calculations since it has a time ordered product 
of currents to which standard perturbation theory applies (T^ u is the amplitude for forward elastic 
photon-nucleon scattering ) 

W^(p,q) = —ImT^faq), T^faq) = i J d A ze iqz (p, nucl|T (J^J^O)) | nucl,p). (2) 



2 For reviews see Refs. 0. 
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The operator product expansion in terms of local operators for a time ordered product of the two 
electromagnetic hadronic currents reads 



i / d A ze i<}Z T(J Ul (z)J U2 (0)) 



N,j 



9v\V2 



QfilQfl2^L,N( ,,2 ' ° s ) 



V\\i.\HV2 



Qu2 QfJ,2 



Q 2 



9v2ii2 ( lv\ < lii\ + 9v\V2 < liL\ ( lv2) ^2,jv( 2 ' a ^ 



^ 3 ...^ JV O J '^ 1 '-'' 1Ar} (0)+ higher twists, 



j=a,V>,G (3) 



where everything is assumed to be renormalized (with /i being the renormalization scale). The use 
of the OPE in the short distance regime ( z — > ) differs from its use in the light cone region. In 
the former case the sum over spin-N extends to a finite value for a given approximation, while in 
the latter (the one we have to deal with) the sum over N extends to infinity. The sum over N runs 
over the standard set of the spin-N twist- 2 irreducible (i.e. symmetrical and traceless in the indices 
Ml> ' ' ' jMiv) flavour non-singlet quark operators and the singlet quark and gluon operators: 



rp\ a ^D^ 2 ■■■D^ N} ip, a 



1,2,..., (n} -1) 



(4) 

(5) 
(6) 



Here and in the following we denote the generators of the flavour group SU{nj) by A Q , and the 
covariant derivative by D^; in addition, it is understood that the symmetrical and traceless part 
is taken with respect to the indices in curly brackets. The functions C J k N (Q 2 / fj, 2 ,a s ) are the 
coefficient functions for the above operators. Since the coefficient functions C k N of non-singlet 



operators depend trivially on the number a (see e.g. Ref. [12] or section 5 of the present article) 
we will use for them the standard notation C k s N . Here and throughout the whole paper we use the 
notation 
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16vr 2 



a s 
Air 



(7) 



for the QCD strong coupling constant. The direct application of the OPE of Eq. (^) to the Green 
function leads to a formal expansion for in terms of the variable q ■ p/Q 2 = l/(2x) i.e. an 



expansion for unphysical x — > oo, 



d 4 ze iqz (p,nucl\T{J^z)J u (0)) | nuclp) 



N,3 



9iiu 




' Ax 2 

9^u + PfiPu-^- 
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2x\ 



Q 2 



A iac\ N( m l/V 2 )+ higher twists, 



where the spin averaged matrix elements are defined as 

(p,nucl|0- ? '' {w '-'^ } |nucl,p) = p {m -P^ n} A{ uc1 N {m 2 Jfi 2 ) 



]=ns,^,G (8) 



(9) 



and m n is the nucleon mass. 

To perform the proper analytic continuation of the representation Eq. (JsJ) to the physical 
region < x < 1 one applies a dispersion relation in the complex x plane to the Green function 
T^y . For electron- nucleon scattering where we have hermitian currents one finds that the Mellin 
moments of the structure functions F k are expressed through the parameters of the operator product 
expansion (||) 



l + (-l) 



N 



M k M = 



l + (-l) 



N ,1 



dx x N - 2 F k {x,Q 2 



i=ns,ip,G 



I 1 



nucl,7V 



(10) 



Please note that the odd Mellin moments of F k are not fixed by this equation. However, all moments 
in the complex N plane are fixed by analytic continuation from the even Mellin moments when all 
the even moments are known. This means that the structure functions in x-space, < x < 1, can 
be found by means of the inverse Mellin transformation when the (infinite set of) even moments 
are known. 

The Q 2 -dependence of the coefficient functions can be studied by the use of the renormalization 
group equations 



2 <} +/?M/x 2 )) 



d 



A* 2 ^ +/?M/x 2 )) 







da s (fi 2 ) 



Cfc,1v ( ~Jj a s (fj 2 ) 



C{ N ^,a s (M 2 )j7^K(^ 2 )) 



i,j=V,G (11) 



Q 



cr N iWM^)) (12) 



where Eq. (|11| ) represents the singlet sector where quark and gluon operators mix under renormal- 
ization, and Eq. (12) is the non-singlet equation. (5{a) is the beta- function that determines the 
renormalization scale dependence of the renormalized coupling constant. It is known at three loops 
[18] in the MS scheme 



da. 



d In// 2 

A) 



(5{a s ) = -(3 a 2 - f3 ia 3 s - 2 a 4 s + 0(a 5 s ), 
(y C A - ^T F n f ) 



4 



Q4 20 
A = (—Cl-AC F T F n f -—C A T F n f ) 

02 = (^C 3 A + 2C 2 F T F n f -^C F C A T F n f 



, >? C\T F n f + -C F T 2 n) + -^CaT 2 ^) (13) 



1415 2 ^ _ , 44^ ^2 , 158^ ^ 2> 
27 

where C F = | and = 3 are the Casimir operators of the fundamental and adjoint representations 
of the colour group SU(3), T F = ^ is the trace normalization of the fundamental representation 
and n/ is the number of (active) quark flavours. The anomalous dimensions 7jv(os) determine the 
renormalization scale dependence of the operators, that is 



din /i : 



:Of'"-" W ' = Of «, i,j=^,G (14) 



d n ns,{m,...,fi, N } _ ~,nsi n \ n ns,{m,...,n N } /irt 

^2^ = -7ivK)^ • (15) 

We define renormalized operators in terms of bare operators as Or = Z Ob and find 

;0 R =(*z)o B = (*z)z-io R => j = -(^z)z-i (16) 



dlnn 2 \dlnfi 2 J \dlnfi 2 J " \dln [i 2 

where it is understood that in the singlet case Z represents a matrix Z lJ . The renormalization 
group equations are solved in the standard form 

CSfrf^,a.(M 2 )) = C fc ">( WQ 2 )) x exp (- (17) 

The solution for the singlet equations has a similar form but since one gets the exponential of a 
matrix of anomalous dimensions one has to define the exponential properly in the singlet case (i.e. 
a T-ordered exponential, see e.g. Ref. [|l9|). Here a s {Q 2 ) = a s (Q 2 / A 2 -^) is the renormalized (i.e. 
running) coupling constant at the renormalization scale Q 2 



Q 2 \ i fh ln ln « 



Ah J ftlnf^ ^ ln 2 ^ 



2 

MS / V MS 



+ 7 — rr ( A 2 In 2 In f ) - (3 2 ln ln f -^1 ) + ft ft - ft ) (18) 



MS 



and A^jg is the fundamental scale of QCD in the MS-scheme. In practice one may use the DGLAP 
evolution equations |2(J for matrix elements of operators at the scale fx 2 = Q 2 (i.e. Q 2 -dependent 
parton distributions) instead of the renormalization group equations for the coefficient functions 
(for perturbative solutions of the DGLAP equations in moment space see e.g. Refs. [||, ^H|). 
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3 The method 



In this section we will discuss the method [22] for the calculation of anomalous dimensions and 
coefficient functions in considerable detail as it applies to the singlet sector. Let us first elaborate 
on some details specific to the dimensional regularization p3| and the minimal subtraction scheme 



[24], and it's standard modification, the MS-scheme |2j, which form a modern basis for multiloop 
calculations in QCD. We use the symbol a s for the renormalized coupling constant and a b for the 
bare coupling constant. Although renormalization constants Z contain poles in e in D= 4 — 2e 
dimensions, anomalous dimensions are finite as D— > 4. This fact gives expressions for the higher 
poles of Z in terms of the first poles of Z. To see this we write Eq. (16) as 

~' z = - (^ z <°- ;>) = " (i z ) 5^? = - ik z ) { - ea - + PMl (19) 

where (3(a s ) is the 4-dimensional beta function of Eq. (|i~3| ) and [— ea + (3(a)] is the beta function 
in 4 — 2e dimensions. This latter function receives no higher order corrections in e due to the form 
of renormalization factors in the minimal subtraction scheme, viz.: 



7 7 i ^° i @ 1 \ 2 , mi 3> 

a b = Z as a s , Z as = 1 a s + (— - — )a s + 0(_ay 

& S — 



d ( a b^ £ ) = = £Za asfI 2e + dZ as da s 2e + ^ da s ^ 2e 
dliifi 2 da din fi 2 3 d In (i 2 

dln/,2 ^a s + Z as 

where a& is a dimensionless object and abfJ- 2e is the bare coupling constant which is invariant under 
the renormalization group transformations. The factors Z %J are calculated as series in a s , and have 
the well known form 

Z v = zv(0) + Z^ l \a s )/e + Z^ 2 \a s )/e 2 + ■■■, i,j = i/>,G 

with = Z GG (°) = 1, Z^ G (°) = Z G ^ = 0. It is helpful to write the matrix equation © as 

4 separate equations, 

ua s 

^ Z QG + ^G Z GG = _ l _ £as + (3(as)] 

oa s 
oa s 

1 Gi > ^G + 1 GG z GG = _ r _ gQs+/?(as)1 JLzGG (21) 
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where we underlined the terms that contribute to the lowest order in e (i.e. order e ). From these 
terms one immediately finds that the anomalous dimensions are expressed through the coefficients 
in front of the first poles of Z 



<} Z^Xas)) i,j = AG (22) 



da 



where Z' L ^ l \a s ) was defined as the order 1/e part of Z l K The coefficients of higher poles in Z can 
then be expressed in terms of 7 iJ by substituting the expression for "f 3 back into equation (^T]). 

The operator product expansion of Eq. (||) is an operator statement and both the coefficient 
functions C l k N and the anomalous dimensions 7^ of the operators are functions and do therefore 
not depend on the hadronic states of the Green function to which one wishes to apply the OPE. 
The information on the hadronic target is contained in the operator matrix elements A l N in Eq. (||) 
which are generally not calculable perturbatively. It is therefore standard to consider simpler Green 
functions with quarks and gluons as external particles, instead of the physical nucleon states, in the 
calculation of coefficient functions and anomalous dimensions. In this case the Green functions can 
be calculated in perturbation theory as well as the operator matrix elements and the anomalous 
dimensions and coefficient functions can be extracted as will be shown below in detail. 

Let us consider the following 4-point Green functions 

T™ J (P, q)=ij d 4 ze^ z (p, quarkjT (J„(*) J u (0)) | quark, p) (23) 

T^{p, q)=ij d 4 ze illz (p, gluon|T (J^z) J„(0)) | gluon, (24) 

where the label 7 is used to indicate an external photon, q indicates an external quark and g an 
external gluon. Spin and colour averaging for the quark and gluon states is assumed. Analogously 
to the decomposition of the hadronic tensor W^ u in terms of F2 and Fl we decompose the Green 
functions in terms of Tl and T2. In the leading twist approximation (i.e. dropping non- leading 
terms in p 2 ) one finds 

n = ww^" T - r2 = -(^i(^^ + ^ 9 i T - (25) 

Applying the OPE to T q7q7 and T g7g7 we find the following equations for the renormalized Green 
functions 



I%™(p,q,a. if J?,e 

y 



00 / -1 \ N 



X 2 •• 2 ' r 



Gk,N( a ^ ~~2 ' e )-^iV^( a s> ~) + G G N (a s , ~T^ £ ) Z N^( a ^ 7) ) ^quark,Ar( a s> 772 ' £ ) 



+ ( Ct^ N (a s , Qj, e)zff G (a s , -) + C G N (a s , Qj, e)Z GG (a s , -) J A G uarkN (a s ,^,e) 

\ Lb o LL c. / Lb 

Tl^(p,q,a s ,^,e) = 



+ 0(p 2 ) (26) 



7 



E 

N=2 



1" V 

2x 



C 



fj, 



[a s ,—,e)ZT • 



1 G Q 2 

a s , -) + C k N (a s , —^,e)Z 



^gluon,7v( a . 



+ lci N (a s ,%,e)Zp(a s ,-J + C? 



1. 



( Q 2 a\V G ^( 



^gluon.JV^S' ..2- 
H 1 



+ o( P 2 ) 



(27) 



where k =2,L, a s = a s (/i 2 /A 2 ) and it is understood that the l.h.s. is renormalized by substituting 
the bare coupling constant in terms of the renormalized one, 



a b 



A) 2 , 
— a s + 



P 2 



2s 



\)al + 0{at) 



(28) 



The terms 0(p 2 ) in the r.h.s. of Eqs. ( ^6| ) and (p^) indicate higher twist contributions. The 
renomalization factors for the external quark and gluon lines are overall factors on both sides of 
the equations and are omited. The coefficient functions on the r.h.s are renormalized quantities. 
The matrix elements A l N are the matrix elements of bare operators and are defined as in Eq. (||) 
with the nucleon states replaced by the appropriate quark or gluon states. 



Figure 1. A graphical representation of Eqs. (|26| ) and (^7j). The symbol <g> indicates 
an appropriate quark or gluon operator (defined in Eqs. (Q) ,(||), @). On the l.h.s. 
of the equations also the crossed diagrams contribute but they are not explicitly 
shown. 
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It is known that the gauge invariant operators and O g mix under renormalization with 
unphysical operators (that are BRST variations of some operators or that vanish by the equations 
of motion) Jl7| ^, But physical matrix elements (i.e. on-shell matrix elements with physical 
polarizations) of such unphysical operators vanish. Since the method that is described below deals 
with physical matrix elements we omited the unphysical operators in Eqs. (^) and (p7|) . 

Starting from Eqs. (p3) and (pTl), the anomalous dimensions and the coefficient functions are 



calculated using the method of projections of Ref. [22]. It reduces the calculation of (moments of) 
coefficient functions and anomalous dimensions to the calculation of diagrams of the propagator 
type instead of the 4-point diagrams that contribute to T^ u . This method relies heavily on the 
use of dimensional regularization and the minimal subtraction scheme and implicitly involves a 
considerable rearrangement of infrared and ultraviolet divergences. 

The method consists of applying the following projection operator to both sides of Eqs. (p6|) 
and (|7]). 

^{mi . . . q^n} qN 



Vn 



N\ <9p« ■ ■ ■ dpw 



(29) 

p=0 



Here q^ 1 ■ ■ ■ q 11 ^} [ s the harmonic (i.e. symmetrical and traceless) part of the tensor q^ 1 ■ ■ ■ q^ N (see 
next section). The operator Vn is applied to the integrands of all Feynman diagrams (nullifying 
p before taking the limit e — > 0, to dimensionally regularize the infrared divergences as p — > 
for individual diagrams). It is important to realize that this operation does not act on the 
renormalization constants Z N and the coefficient functions on the r.h.s. of Eqs. (p6|), (|27j). It does 
however act on the matrix elements A N . The nullification of p has the effect that of all the diagrams 
that contribute to the perturbative expansion of A l N only the tree level terms (i.e. with no loops) 
survive since all diagrams containing loops become massless tadpole diagrams. Massless tadpole 
diagrams are put to zero in dimensional regularization. Furthermore, the iV th order differentiation 
in the operator Vn has the effect that Vn projects out only the N th moment since of all the factors 
l/(2x) N ' only l/{2x) N gives a non zero contribution after nullifying p. On the left hand side the 
effect of Vn is to effectively reduce the 4-point diagrams that contribute to to 2-point diagrams 
(this follows from the nullification of the momentum p), which drastically simplifies the calculation. 
We apply the operator Vn after the tensor structures 2, L have been projected out because the 
operator Vn would mutilate the tensor structure of T^ u . In the projector Vn we use the harmonic 
tensor q^ 1 ■ ■ ■ q^ N ^ to remove higher twist contributions (the 0(p 2 ) terms in Eqs. ( p^ ) and (|2j 
that after differentiation with respect to p^ survive as terms proportional to the metric tensor. 
Summarizing, we have after application of the projection operator Vn to Eqs. (p^) and (E~ 



n^(^,a s ,s) = (cl N (a s , ^,e)Z^{a s , \) + C% N (a s , ^,e)Z^(a s , i)j A^ N (e) (30) 



T^(%,a s ,e) = (ct N (a s ,%,e)zt f G (a s ^) + C G N (a s ,%,e)Z GG (a s , V) A^N^) (31) 
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where k = 2, L and we denned 



Q 2 q {^...g^} QN 

T k;N (—,a s ,s) = - QpfMi dpm T k (p, q ,a s ,v ,e) 



(32) 

p=0 



It should be understood that (^) and represent a large coupled system of equations when 
both sides are expanded in powers of a s and e (i.e. C is expanded in positive powers of e and Z is 
expanded in negative powers of e). 

After the calculation of T q7q7 and T g7g7 in the order and the determination of the tree 
level matrix elements Ajy ree one can solve Eqs. ( |30[ ) and (|3l]) simultaneously to obtain , C G , 
and Z^ G in order but, unfortunately, Z G ^ and Z GG only in the order a 2 . This limitation 
follows directly from the fact that starts from the order a® but C G starts from order a s since 
the photon couples directly only to quarks. In solving the equations it is essential that all poles 
of the Z factors are fully expressed in terms of the anomalous dimensions as was discussed in the 
beginning of this section. Coefficient functions and operator matrix elements are finite as e — > 
but one must make sure that sufficiently high powers in e are taken into account. For example, 
one should consider order e 2 contributions for Ct at order a s . We stress that by calculating only 
propagator type diagrams in the l.h.s. of Eqs. (|30|) and ( |3l| ) we can get both renormalization 
constants of operators and coefficient functions. 



Figure 2. Examples for diagrams contributing to the Green functions T g ^ g< ^ (a) and 
T<i<t><i<f> (b). 



To obtain Z G ^ and Z GG in order we calculated two more unphysical Green functions T q<?iq ^ 
and r g ^ s ^ (see fig 2), in which the photon is replaced by an external scalar particle (f> that couples 
directly only to gluons. The vertices that describe the coupling between the external scalar field 4> 
and the gluons follow from adding the simplest gauge invariant interaction term 4>G'^ iy G% u (where 
G®„ is the QCD field strength tensor) to the QCD Lagrangian. For the Green functions r q< ^ q<?i and 
Ts4>s<P an OPE similar to @ exists with the same operators but with different coefficient functions 
C G and C^, where starts from the order a s and C G starts from the order 

Repeating the steps that led to Eqs. ( |30| ) and ( pl| ) one finds for these Green functions the 
following equations 

{Z G ,) 2 Tf^{%,a s: e) = (c G N (a s , %,e)Z GG (a s , \) + ^ N (a s , %,e)Zf(a s , h) A^ >N (e) 

LL \ LL c Lb c. / 

(33) 
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LL \ LL c. LL o / 

(34) 

As is indicated in the l.h.s., the external operators Gp U G% u have to be renormalized [26| 



[G^ G aJ R = Z G^ B H > Z G 2 



1 



fl M'-Vb 1 ' 1 - /3(a.)/(a 8 e) 

where the dots indicate (unphysical) operators that are omited since they do not contribute to 
the on-shell matrix elements with physical spin projections that we consider. (The only physical 
operator that mixes under renormalization with G^G^f is the quark operator uiqipip that vanishes 
in our limit of massless quarks.) We emphasize that the 3- loop (3- function is required for the present 
3-loop calculation. 

After the calculation of T q< ^ and T s< ^ in the order a\ one can solve Eqs. flU) and @ to 



obtain C G , Cj, Z G ^ and Z GG in the order a 3 s (from these equations Z^ G and can be obtained 

in the order a 2 only). Please note that the coefficients C G and C"f are obtained as a byproduct 
and are not important for the physical process under consideration. Furthermore, the two sets of 
equations (|0|), (gj) and (||), @ deter mine all the anomalous dimensions of the order a 2 in two 
independent ways which provides a consistency check on the results. 

Solving the system of equations for Fl is slightly more involved than for Fi since the structure 
function Fl contains no tree level (i.e. order a®) contribution. To solve the sets of equations for 
Fl one should add extra information, for example the 0(a 2 ) contributions to Z G ^ and Z GG as 
determined from the equations for F2. 



4 The calculation 

As was discussed in the previous section, we will apply the operator Vn for ^=2,4,6,8 to 4 different 
Green functions, T^ 7q7 , T^ 7g7 , T q ^ q ^ and T g ^ g< ^ and we sum over the physical spin polarizations 
of the external quarks and gluons. For the external quarks (in T q7q7 and T q?iq ^) the sum over 
the polarizations is performed by inserting the projection operator between the external quark 
legs and taking the trace over the strings of gamma matrices. For the external gluons (in T^ 7S7 
and T g< ^ g< ^) the sum over physical spins can be done by contracting the external gluon lines with 
—g a/3 + (p a q l3 +p l3 q a )/p-q—p a pPq 2 /(p-q) 2 in which the (on-shell) gluon has momentum^ (withp 2 = 
0). The presence of the extra powers of p poses considerable efficiency problems (the operator Vn 
will generate more than 3 times larger intermediate expressions as compared to the case of a simpler 
gap projection). Alternatively one may take the sum over physical gluon spins by contracting the 
external gluon lines with only —g a ^ and adding external ghost contributions to the Green functions, 
T^ 7h7 to Tf 7g7 and T h< ^ to T g ^, where the label h indicates an external ghost line. This 
procedure is identical to the standard use of ghost diagrams to remove unphysical polarizations 
of gluon propagators in the covariant gauge, and the ghost particle h is the same ghost that we 
use in closed loops. Although we now have to consider all diagrams that contribute to T^ h ^ and 
rph(f>h<p^ ano i i ncrease total number of diagrams that we have to calculate, it still makes the 
computations more than a factor of 3 faster (since ghost diagrams are of a far simpler nature than 
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gluon diagrams). We checked for the lowest moments that the two methods for taking the sum over 
the gluon spin polarizations gave the same results, but for the higher moments we only applied the 
ghost method. 



Figure 3. Examples for ghost diagrams contributing to the Green functions T h 7 h 7 ( a ) 
and T h< ^ (b) 



The explicit generation of Feynman diagrams with the corresponding symmetry factors has 



been done automatically by the use of the program QGRAF[27]. Statistics on the number of 
diagrams in the different classes q7q7, q<M</> etc is presented in table 1. The generation (and the 
counting) of the diagrams is non-standard, since a number of tricks have been used (for example, 
crossed diagrams are not generated explicitly, and the diagrams with 4-gluon vertices are split into 
a number of parts for which the colour factor necessarily factorizes) . 





Tree 


1-loop 


2-loops 


3-loops 


Lorentz projections 




1 


3 


27 


413 


2 






1 


24 


697 


1 


g7g7 




2 


20 


366 


2 


117117 






2 


53 


2 


g4>g(j) 


1 


11 


241 


7219 


1 


h(f)h(f) 




1 


36 


1266 


1 


TOTAL 


3 


23 


399 


10846 





Table 1. Number of diagrams and Lorentz tensor structures in the classes q7q7, 
q4>q4>, g7g7, hyhy, g(j)g(j) and hc/dup. Notation: q = quark, g = gluon, h = ghost, 7 
= photon, <j) = scalar particle that couples only to gluons. 

It is clear from these statistics that the calculation of the diagrams necessarily has to be automated 
to a large extent. The calculation is therefore organized as follows: 

1. The diagrams are generated automatically with a special version of the diagram generator 
QGRAF[^?]]. For every class q7q7, q</>q^ etc. the full set of diagrams is put into a single file 
using a dedicated database program MINOS that manages information about thousands of 
diagrams and can be instructed to call other programs, giving them the proper information 
from a database file. 
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2. The representation for a diagram at this point is still a very compact one, and this is explored 
as follows. We use programs written for the symbolic manipulation program FORM to 
calculate colour factors for each diagram and bring the diagrams into a representation that 
explicitly contains information required at later steps in the calculation. For instance, this 
involves choosing automatically an optimal path (in most cases the shortest path) for the 
external momentum p to flow through each diagram (we are going to expand in p when the 
operator Vn is applied) and determining automatically the diagram's topology when p is 
nullified. This information, for all diagrams together, is kept in a single file and is accessible 
using MINOS. 

3. We instruct MINOS to run sequentially, one diagram at a time, a highly optimized FORM 
program that performs the explicit calculation i.e. it substitutes all the Feynman rules, it 
performs projections on the Lorentz structures of the Green functions, it Taylor-expands the 
diagram in the external momentum p (the depth of the expansion increases with the moment 
index N), it takes all the Dirac traces, contracts with the tensors q^ 1 ■ ■ ■ q^ N ^ and finally calls 
the MINCER |^] integration package to perform the 3-loop scalar integrals of the massless 
propagator type (using the integration by parts algorithms published in Ref. pC|]). The 
results together with some useful technical information about the calculation (such as the 
resources used) are again stored into a single file. MINOS will initiate the calculation of a 
next diagram as soon as the calculation of a previous diagram is completed without any need 
for human interference. 

An important aspect of the calculation is the use of the symmetrical and traceless tensors, 
qivi . . . gMiv} that are used to extract the leading twist contributions. These tensors are (with the 
proper normalization) known as 'harmonic tensors' H n satisfying 

H^-^S^ 2 = (traceless) 
//,'; : "''' = /AT"'' (symmetrical) 

H^-^QtMn = H^^Q 2 (normalization) 

and an explicit construction (in Euclidean space-time) in terms of the Kronecker delta symbols S^" 
and the momenta reads 

h^= h l(Q 2 y /2 E 8(j*i,-~>f*j)Q lii+1 ~-Q ,tn , 

j=0,2, --,n index perm. 

h n_ ( /2 /2 T(2-2e)T(l-e + n-j/2) 

h * ~ 2 r(2-2e + n)r(l-e) ^ 

where the second summation is over all the ways to partition the indices into one set containing 
j indices (put into 5 ) and a second set containing (n — j) indices (put on momenta Q ■ ■ ■ Q). T 
is the Euler gamma function (factorial function). The tensor 5 is a completely symmetrical tensor 
constructed from Kronecker delta symbols only, 

<$(/ii,/i 2 ,/^4) = S^S^ + S^S^^ + d^^S^ 3 , 
etc. 
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such that the normalization is 

u n on r ( 2 " e + n ) 

0(/il,Ml,/i2,At2 ; • • • ,Mn,/V) = 2 _ g ^ 

For example the tensor with four indices is 



-hlQ 2 



gHalUQIHQM + §t*2t*4Q^Qt*3 + §1*21*3 QWQW 



A' 4 



+ ^1M4QM2QM3 + j/iiWQWQW + §t*it*2QI*3Q 

The number of terms in a harmonic tensor increases rapidly with the rank of the tensor and in 
the present calculation these harmonic tensors are contracted with many other tensors constructed 
from only a few different (integration) momenta. An efficient implementation of the harmonic 
tensors that directly takes into account the symmetries of the other tensors can greatly limit the 
number of terms that are produced in the calculation and is vital for its feasibility^. 

The results of our calculations are presented in the last section of this article. The results for 
the second moment N=2, that are presented in this paper, were obtained in an arbitrary covariant 
gauge for the gluon fields. This means that we keep the gauge parameter £ that appears in the gluon 
propagator i [—g^ + (1 — CjQ.^Q 1 ' /(l 2 + 2e)]/((7 2 + ze ) as a f ree parameter in the calculations. The 
explicit cancellation of the gauge dependence in the coefficient functions and anomalous dimensions 
gives an important check of the results. 

Another non-trivial check of the calculation is the validity of the singlet relations 72 + 12 ^ = 
0, 72"^ + 72^ = that follow from the fact that the QCD energy momentum tensor (which 
contains both the spin-2 quark and gluon operators) is conserved and has therefore zero anomalous 
dimension. 

Please note that at present the only independent check of our 3-loop results is provided by the 



calculation in Ref.[31] where the leading nj terms were calculated for the non-singlet anomalous 
dimensions in all orders in a s using a large n/ expansion. Our afn^ terms agree with Ref.Q. 

The reader may notice the appearance of the constant £4 in the final 3-loop results for the 
coefficient functions C^^ and which seems to be in conflict with the empirical law that in the 
results for inclusive physical quantities £4 does not appear. However these 3-loop coefficient func- 
tions contribute to the next-to-NNL (NNNL) order for F2 and require 4-loop anomalous dimensions 
to get a complete physical NNNL approximation. Assuming the cancellation of £4 in the complete 
NNNL approximation, one can derive the coefficients of £4 in the 4-loop anomalous dimensions. 

In spite of all optimizations of the integration program MINCER and of many other efficiency- 
crucial parts of the calculation, the calculation of the singlet moments, A^=2,4,6,8 as published 
in this paper, required more than the equivalent of 7500 hours on a 150 Mhz SGI Challenge 
workstation and required at some instances 2 Gbyte of storage place for the intermediate stages 



3 it is therefore interesting to mention that an efficient implementation of Eq. (|3fj) exists in FORM-2 and requires 
only 3 lines of code using the built in combinatorical functions. 
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in the calculation of one of the diagrams^]. We noticed that both the required disk storage space 
and computation time increased with almost a factor 5 when we compared the calculation of each 
(N+2)th moment to the calculation of the corresponding Nth moment. The calculation of higher 
moments using the same methods is therefore not feasible at present. 



5 Analysis 

In this section we will investigate the effects due to the calculated 3-loop coefficient functions and 
anomalous dimensions (the list of the analytical 3-loop results is given in the next section). We 
will reconstruct the coefficient functions and splitting functions as distributions in x-space using 
the calculated moments and incorporating as much as possible the known information about the 
leading (and in some cases the next-to- leading) singularities in the limits x — > 1 and x — > 0. Explicit 
comparison with the known 2-loop coefficient functions and anomalous dimensions shows that our 
approach to the reconstruction gives very good effective fits of the x-space distributions when only 
few moments are used and when sufficient information about the two endpoints x = 0, x = 1 is 
known. This procedure also allows us to estimate the error of the reconstruction. 

In order to relate the calculated results for the anomalous dimensions 7jv and the coefficient 
functions Cat to experiment one must be able to obtain the experimentally measurable structure 
functions F^{x, Q 2 ) from the Mellin moments M^j^ of Eq. (|To|). The rigorous procedure to obtain 
Fk(x, Q 2 ) from the moments is to apply the inverse Mellin transform which, however, requires the 
exact knowledge of jy in the complex iV-plane (or equivalently the analytic continuation to 
complex N from all even or odd moments). Because at the 3-loop order we have calculated only 
a limited number of moments jn and Cat, we can obtain only approximate results in x-space. An 
example of a NNL order analysis in x-space based on a limited number of non-singlet moments can 



be found in Ref. |3J. 



As an alternative to studying the Q 2 -dependence of (x,Q 2 ) via the Q 2 -dependence of moments 
Mk,N (see Eq. (O)), one may start from the x-space distributions j(x) and C(x) that are related 
to 7iv and Cn via the Mellin transform 



IN 



dxx N - 1 -f(x), C N = dxx N ~ 1 C{x) (36) 



and do the Q 2 -evolution in x-space via the DGLAP equations 



2 d ( q+(x,Q 2 ) \ _ ( i++(x,a.) ^ G (x,a s ) \ ( q+(x,Q 2 ) . 




Ql ' ' ! l ns > + (*, a.) U{x, Q 2 )- q -^l^l) (38) 



(qj(z, Q 2 )) = - 7 ns '"(x, a s ) ® qj(x, Q 2 ) (39) 



4 This is called 'intermediate expression swell'. It is a well known phenomenon in Computer Algebra. 
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where = qf ± qf, q + = 2/=i 9/ i Qf £ {n, d, s, c, t, 6} are the quark distributions and G is the 
gluon distribution. It is understood that the inverse Mellin transforms of the anomalous dimensions 



are related to the standard Altarelli-Parisi splitting functions P l3 (x) as ^(x) 
symbol (g> indicates the convolution integral 



.put 



f 1 f 1 f 1 dx^ x 

A(x)®B(x) = / dxi / dx 2 A(xi)B(x 2 )5(x - xix 2 ) = / — -A( — )B{x x ). 

JO JO Jx Xi X\ 



The 



(40) 



Since we calculated the moments only for electromagnetic scattering we are restricted to the "+" 
type anomalous dimensions and coefficient functions for the non-singlet sector, i.e. to the distri- 
butions that correspond to analytic continuation from the even moments (i.e. we do not consider 
the "-" type distributions which correspond to analytic continuation from the odd moments and 
which are relevant for processes that are not symmetric with respect to crossing symmetry such as 
neutrino scattering via VF-boson exchange). The proper relation between the structure functions 
and parton distributions for electromagnetic scattering is 



F k (x,Q 2 ) 




C™> + (x,a s )®(af(x,Q 2 )-±-q + (x,Q 2 ) 



1 



C%(x, a 8 ) <8> q + (x, Q z ) + Cg(x, a a ) ® G(x, Q z ) 



(41) 



where et is the electromagnetic charge of a quark of the flavour /. 

A special feature of the Q 2 -evolution in x-space (that is a direct consequence of the convolution 
integral) is that the structure functions F k (x, Q 2 ) at some value x = xq only depend on the parton 
distributions, splitting functions and coefficient functions in the region xq < x < 1. This means 
that one does not need to know the anomalous dimensions and coefficient functions below the 
experimentally accessible x-region (but one always needs them for x ~ 1). 

One can see that, away from the singular limits x — > and x — > 1 that received special attention 
in the literature, the 1-loop and 2-loop j(x) and C(x) Q-jll] behave smoothly (as a smooth 
interpolation between the small-x and large-x regions). As a consequence the distributions j(x) 
and C{x) can be approximated to a high precision by linear combinations of simple distributions 
that contain the singular^ terms of the expansions of the exact j(x) and C(x) in the limits x — > 
and x —s- 1, i.e. log l (x)/x and log l (x) for small x and [log l (l— x)/(l— x)]+, log l (l— x) and 5(1— x) for 
large x plus finite order polynomials of x (the order determines the accuracy of the approximation) . 
(The exact results for 7(2;) and C(x) contain in general complicated polylogarithmic functions but 
after their expansions in the small-x and large-rr regions one finds the simple structure of the 
singularities mentioned above.) 

The most singular terms of j(x) and C(x) in the limit x — > (i.e. log l (x)/x ) show up as 
singularities in moment space for N = 1. By increasing the moment index N, one increasingly 

5 [/(a;)]+ type distributions appear below, where fix) is non-integrable at x = 1. The standard definition of 
integrals over these distributions is: J* g(x)[f(x)]+dx — (g(x) — g(l)) f(x)dx — g(x)f(x)dx. Or in a more 
symbolic form [f(x)]+ = f(x) - 5(1 - x) J Q f(z)dz. 
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probes the region x ~ 1 (since the factor x 1 suppresses the contribution from the small- x 
region). In particular, when the moments grow for large N this is due to terms of the type [log*(l — 
x)/(l — x)]+, % > (of which the moments increase as log* +1 (iV) for large N). 

The approximations of 7(x) and C{x) in terms of the linear combinations of simple distributions 
(i.e. log l (x)/x, log*(x), [log l (l — x)/(l — x)] + , log l (l — x), 5(1 — x) and low powers of x ) provides 
a natural scheme to study the effects of the 3-loop moments N=2,4,6,8,(10) on the Q 2 -dependence 
of structure functions in x-space and to quantify the uncertainty due to our lack of knowledge of 
the higher (i.e. N > 10) 3-loop moments. Considering such approximate distributions for the a 
priori unknown j(x) and C(x) at 3-loops we may impose the correct asymptotic moment behavior 
for N — > oo and N — > 1 by including in the linear combinations of simple distributions the correct 
leading (and subleading) terms in the large- x and small- x regions. In some cases information about 
these leading (and subleading) terms in the large and small x regions is available in the literature. 

The linear combinations of distributions should also include various functions of a less singular 
type (compared to the leading singularities), e.g. log l (x), log l (l — x) and low powers of x, which are 
relevant terms in the intermediate x region. We fit all unknown coefficients in a linear combination 
of simple distributions such that the moments N=2,4,6,8,10 of the linear combinations give the 
calculated values. Since we know only 4 singlet (and 5 non-singlet) moments explicitly we can 
only allow up to 4 (or 5) arbitrary coefficients (i.e. types of simple distributions) in the linear 
combinations. Only in the case of the quark anomalous dimension we can allow more than 5 
arbitrary coefficients by considering pure-singlet (ps) and non-singlet (ns) parts separately (since 
information about the small-x and large- x behaviour is known for these parts separately). 

Although at the 3-loop order there are only a finite number of allowed simple distributions of the 
types log l (x)/x, [log l (l — x)/(l — x)] + , 5(1 — x), log*(x), log l (l — x) we do not have enough moments 
to include all of them simultaneously in a fit and we have to choose among the possible simple 
distributions. Therefore, to obtain an estimate of the uncertainty of the reconstructions due to our 
lack of knowledge of the higher moments we studied the stability of the x-space reconstructions 
with respect to a change of the simple distributions included in the linear combinations. If one 
does not find stable results under such a change one should have included a larger number of simple 
distributions of the types mentioned above. However, as one can see below, in many cases we have 
sufficient information to obtain stable results. 

For a correct x-space reconstruction it is crucial that one knows the moments exactly (or with 
a high precision) since the inversion of moments into x-space distributions from a limited set of 
moments is numerically unstable, i.e. a small change in one of the moments can lead to a large 
change in the reconstructed result. In relation to this, we note that also the dedicated procedure 
for reconstructing 7(2;) and C(x) that we consider here requires very accurate numerical values 
for the moments. Therefore one can not rely on interpolation between known even moments to 
increase the number of available moments. 

5.1 Anomalous Dimensions 

We will now consider the 3-loop singlet and non-singlet anomalous dimensions. We will use the 
following notation for the inverse Mellin transformed anomalous dimensions 

7 ns (x,a s ) = a sl ns ' {0 \x)+ah ns ' W (x)+a 3 s Y ls '^(x) + O(a 4 s ) (42) 
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f j (x,a s ) = a s f^ Hx)+ah tUl \x)+ah tm (x)+O(aj), i,j = ^,G (43) 



where a s = a s /(47r). In Ref, [|33| it was proven that non-singlet anomalous dimensions are finite 
as N — > 1 which means that in x-space the non-singlet anomalous dimensions are less singular 
than 1/x as x — > 0. The leading singularity has been derived pi] from the leading order small- x 



resummation of the non-singlet evolution kernels [35] 



7 W2) (x ^ 0) = KjZ log 4 (x) + (\og\x)). (44) 
The singlet anomalous dimensions are in general singular as N —* 1 ( ie x — > ). The following 



results, which have been derived using small- x resummations, may be found in the literature [36 
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yW>,(2) (x _ 


-0) 


7 VG,(2) (x _ 


-o) 


7 G ^' (2) (x - 


-o) 


7 GG,(2) (x _ 


-o) 



896 

: — n f C A C F \og(x)/x + 0(l/x) 

= ^ nf C A log(x)/x + 0(l/x) 
- 0(\og(x)/x) 

■■ O (log(x)/x) . (45) 

The last two equations simply mean that there are no leading singularities of the type log 2 (x)/x 
(since the term in the BFKL anomalous dimension [|36| vanishes). 

For the limit x — > 1 of the anomalous dimensions not much is known except for the conjecture 
0] that the diagonal elements 7^ and 7 GG do not rise faster than log(iV) as N — > 00, which means 
that one has in x-space a leading term [1/(1 — x)] + (and not [log(l — x)/(l — x)] + ). This conjecture 
is based on the explicitly known one and two loop results. But the rise of the 3- loop non-singlet 
moments observed for N = 2,4, 6, 8, 10 indicates that it also holds at the 3-loop level (putting both 
the terms [1/(1 — x)]+ and [log(l — x)/(l — x)]+ in the linear combination of simple distributions 
gives a small numerical coefficient for the second term compared to the first one, although this is 
clearly not a proof). Furthermore, the one and two loop coefficients of the terms [1/(1 — x)]+ in 
7^ and 7 GG have a simple ratio Cf/Ca and we presume that this ratio is the same at 3-loops. 

The off-diagonal singlet anomalous dimensions 7^ and to zero as N — > 00. The best 

way to see that they do not contain terms [log*(l — x)/(l — x)]+, i > 0, is to consider the approach 
of Refs. [Q, IH, |(| to calculate the anomalous dimensions, i.e. to study the renormalization of the 
singlet twist-2 operators where one can see that the terms [log*(l — x)/(l — x)]+ originate from 
diagrams that only appear for diagonal anomalous dimensions (notice that this method allows the 
direct calculation of these terms without performing the complete calculation of the anomalous 
dimensions). Therefore, the leading terms in 7^ and 7 G ^ as x — > 1 that one can expect from 
performing the necessary integrals at 3-loops are log 4 (l — x). 

Finally we presume that the pure singlet (ps) combination j ps = 7^ — j ns rapidly vanishes 
as N — > 00 (i.e. vanishes at least as quick as 1/N for large N; at 2-loops it is known to van- 
ish as 1/N ). We have observed this tendency already in the low-N moments also at 3-loops. 
This means that 7^ and r ) ns contain the same terms that are important for x ~ 1 (since the 
transform from x- to N-space for large N gives [log*(l — x)/(l — x)] + — > log l+1 (N), 5(1 — x) — > 1 
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and log*(l — x) — > log 1 (N)/N ) and that they only differ in terms that are important for small- x 
In relation to this, recall that only 7^ contains a term log(x)/x and j ns does not, see Eqs. (fS|,45) 



Figure 4. The exact 2-loop anomalous dimensions and reconstructions (approxs.) 
based on the lowest 5 moments N=2,4,6,8,10 (rt/ = 4). The approximations are 
obtained by fitting the sets of distributions of Eq. (|46|) to these 5 moments 
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At 2-loops the analysis based on only a few low-N anomalous dimensions can be compared with 
the exact results as is shown in figure 4. The indicated approximations in figure 4 are based on 
matching the linear combinations of simple distributions of Eq. ( |4"6| ) to 5 moments (N=2,4,6,8,10) 
of each anomalous dimension j ps , j ns , 7^ G , 7 G ^ and j GG . 



7 V>G,(1) )7 GV>,(1) 



7 



GG,(1) 



7 



ns,(l) 



1/x, log(x), 1, log(l — x), log 2 (l — a?)} 
l/x,log 2 (x), l,x,log 2 (l - x)| . 

l/x,log 2 (x),l,log 2 (l -x 
l/x,log(x),l,log(l -x), 
log 2 (x),log(x),log(l - x) 



1-x 



,6(1 -x) , 



,6(1 



1-x 



,6(1 



log 2 (x),l,log 2 (l-x) 

1/x, log 2 (x), log(x), 1, x j 
1/x, log 2 (x), 1, x, x 2 | 



,6(1 



( 4 ?) 

In these formulae each pair of curly brackets encloses a set of simple distributions that appear in 
a linear combination with coefficients to be matched. We take 2 different linear combinations for 
each anomalous dimension to see the stability of the fit. 

Please note that all the coefficients in these linear combinations were taken as arbitrary in the 
fit (i.e. we even did not fix the coefficients of the leading singularities) to have a strong check of the 
approach. The only additional information that is used in the reconstruction is that the coefficients 
of the most singular large- a; terms [1/(1 — x)]+ in 7^ and 7^ have a simple ratio Cp/Ca- 

The considered distributions are singular at x = and x = 1, and because e.g. 6(1 — x) can 
not be drawn one should consider the figures only as an indication of the accuracy of the full 
distributions which are used in convolutions (where all the terms that are singular at x = 1 fully 
contribute). The approximations are stable for x > .2 and do dramatically improve in the small- 2; 
region when more moments (N=12,14,- • •) are used such that more simple distributions can be 
included in the linear combinations. 
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Figure 5. The effects of the NNL order (i.e of the 3-loop anomalous dimensions and of the 
3-loop beta-function) on the evolution of the singlet distributions q + (x,Q 2 ) and G(x,Q 2 ). The 
parton distributions q^ L , Gnl and q^ NL , Gnnl are obtained as solutions of the DGLAP evolution 
equations ( |37| ) in respectively the NL and NNL approximations starting from the parametrization 
scale Qq = 4 GeV 2 . The difference between curves for the same Q 2 indicates the uncertainty in the 
NNL order effects due to the lack of knowledge of the higher anomalous dimensions 7V, N > 8. 
The parton distributions at Qq correspond to the MRS (A) set [38|. For simplicity, nj = 4 and 



Aqcd = 300 MeV are taken throughout the evolution. 

The effects of the NNL order on the evolution of the singlet distributions q + (x, Q 2 ) and G(x,Q 2 ) 
is illustrated in figure 5. The evolutions are done using a program based on the Laguerre polyno- 
mial technique of Ref. []39|]. An estimate of the uncertainty in the NNL effects due to unknown 
higher moments is made by fitting various linear combinations of simple distributions to the cal- 
culated 3-loop moments of the anomalous dimensions. The curves correspond to the sets of simple 
distributions of Eq. ( |47| ) (two sets for each anomalous dimension) which are representative for a 
larger variation in the choice of simple distributions. 



7 V>G,(2) 



7 G<M2) 



log(x)/x,log(x),l,log(l - x),log 4 (l - x)| , 
log(x)/x, 1/x, log(x), log(l — x), log 3 (l — x)\ 
log(x)/x,log 2 (x),log 2 (l - x),log 4 (l - x)| , 
\log(l - x),log 2 (l - x),log 3 (l - x)} , 
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7 



GG,(2) 



l,log 2 (l-x),log 4 (l-x), 



log(x)/x, l,log 2 (l - x), 



1-x 



1-x 



+ ,8(1- 
,8(1 -x) 



7 



ns,(2) 



1-x 



7 



log 4 (x),l,log 2 (l-x),log 4 (l-x), 

log 4 (x),log 2 (x),l,log 2 (l -x), 

,log(x)/x, 1/x, log 3 (x), 1, x 2 | 
[ log(x)/x, log 4 (x), log 2 (x), 1, x 2 } 



1-x 



,8(1 



,8(1 



The coefficients of the leading x — > terms of ^ G '( 2 ), ^p s >( 2 ) and j ns '^ are known exactly, see Eqs. 
(^ , ^5|) . Please notice that in the variation of simple distributions we changed even the type of 
the expected leading singularities as given in Eq. ( f47| ) (when little information about these leading 
singularities is known) to obtain a conservative estimate of the accuracy of the reconstructions. 
One may see that the NNL effects are of the order of a few percent in the region of x > 0.1 . 
Approximately 50% of the shown NNL order effect on the evolution of the parton distributions is 
due to the inclusion of the 3-loop beta-function and the other NNL order terms in the expression for 
the strong coupling constant, Eq. (|l8|). The other 50% is due to the 3-loop anomalous dimensions. 
The estimated uncertainty is small in the evolution of q + (x) and is somewhat larger in the evolution 
of G(x). This is a consequence of the especially good reconstruction of the quark anomalous 
dimensions (see figure 6), which is itself due to the existence of separate information about the 
pure-singlet and non-singlet parts. 

The curves in figure 6 correspond to matching the following linear combinations of simple 
distributions to the calculated moments (where the coefficients of the leading x — > terms are 



known exactly, see Eqs. (44j45|)) 



7 



ns,(2) 



log 4 (x),log 2 (x),l,log 2 (l-x) 



1-x 



7 



log 4 (x), log(x), log(l — x),log 3 (l — x) 

log 4 (^),l,log 2 (l-x),log 4 (l-x), 

log(x)/x, 1/x, log 4 (x), log 2 (x), 1 j 
log(x)/x, 1/x, log 3 (x), 1, x 2 j 
log(x)/x, log 4 (x), log 2 (x), 1, x 2 j 



1-x 



,8(1 -x)\, 

8(1 -x] 
8(1 -x)\, 



1-x 



(48) 
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Figure 6. The quality of the reconstruction of the 3- loop quark anomalous dimensions based on 
the moments N=2,4,6,8,(10 for non-singlet only) and the known leading small- x terms for -yP^ 2 ) 
and 'y ns '( 2 \ The reconstructions are obtained by fitting the sets of distributions of Eq. (|4g| ) to the 
available moments. 

From the foregoing analysis we have concluded that the inclusion of the 3-loop anomalous 
dimensions in the evolution of parton distributions in the MS-scheme changes both the quark and 
gluon distributions by 1-3 % in the region 0.1 < x < 0.8 for a change in Q 2 between (2 GeV) 2 and 
(100 GeV) 2 . 

We further remark that it is important to obtain the exact results for the 3-loop anomalous 
dimensions for all N to do a complete analysis in the NNL order. 



5.2 Coefficient Functions 



As an example of the effects of the calculated 3-loop coefficient functions on the deep inelastic 
structure functions we consider the case of CV; S ■ We use the following notation for the inverse 
Mellin transformed coefficient function. 



a.) = C^\x) + a s C™' {l \x) + a 2 s C™' {2 \x) + a^Cf ,{3) (x) + 0(a 4 s 



(49) 



Please note that for the complete NNNL order approximation for F2 one would need, besides the 
3-loop C™'^, the 4-loop anomalous dimensions which are difficult to calculate at present. 

Presuming that the anomalous dimensions are not more singular than log(iV) as N — > 00 one 
has shown that all the logarithms log* (iV) that are present in the non-singlet parts of the coefficient 



functions exponentiate according to the soft gluon resummation formulae [40, 41, 42, 43]. We will 
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use this to obtain the large N limit of moments of the 3-loop CV; s coefficient function (i.e. the x — ► 1 
limit of C™)- Exponentiating the log J (iV) terms from the explicitly known 1 and 2-loop coefficient 
functions we find 



« exp< 



a s [2C F log 2 (N) + C F (3 + 4 7E ) log(iV) 
+O(log 2 (A0)l + <D(a 3 s log 4 (N))\ + 0(a s N° 



2 -C F {^C A - 2 -n f )\og\N) 



(50) 



where ~/e is the Euler constant. In x space this corresponds top 



Cr{x) X ^ 1 \ + a s {AC F 
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3C F 
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-(-fC A C F + \C F n s - 18C F ) 
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+ 0(6(1 -x))\ +a 2 s {8C F 
+ 



log 3 (l-^) 
1-x 



log(l-a) 
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+aUSCl 



3 (i-^) 
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-O 



log 3 (l— x) 
l-x 



+ + (-?fC A C F + §C F n f -30C F ) 
+ o(a 4 



log 4 (l-z) 
l-x 



1-x 



(51) 

The non-singlet coefficient function C ns does not have 1/x terms for small-x (it is finite as 
N — > 1) and the leading small-x term that one may expect at 3-loops (from extrapolating the 
1-loop and 2-loop results) is log 5 (x). (The singlet coefficient function contains the same large- 
x terms since it contains the same non-singlet contribution. The small-x limits of the singlet 
coefficient functions do contain terms of the type \og % (x)/x, i > 0. At 3-loops one expects a term 
\og(x)/x.) 

At the two loop order the reconstruction of the coefficient functions can be compared with the 
exact results. In figure 7 we illustrated the quality of the reconstruction based on matching the 
following sets of simple distributions to the first 5 moments N=2,4,6,8,10 
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r log 3 (i-x)' 
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rio g 3 (i-a-)l 
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+ ' 


l-x 


+ ' 


l-x 



( 52 ) 

It is understood that for the first set the coefficients of the last two terms are taken from Eq. 

was taken exactly. 



fl5lD , and that for the second set only the coefficient of [log (1 — x)/(l 

ns (2) 

One can see from figure 7 that the reconstruction of the distribution C 2 is quite stable for 
x > 0.2 . We want to emphasize that considering only the two most singular soft gluon terms (i.e. 
[log 3 (l — x)/(l — x)] + and [log 2 (l — x)/(l — x)] + ) gives a huge overestimate of the true effects. 

6 for a table of the moments of the relevant distributions in the large TV limit, see e.g. Ref E^]. 
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However, as it was observed in Ref. [1C], all the terms of the type [log l (l — x)/(l — x)]+, i > 0, 
together constitute a large fraction of C^ s A in the region x > 0.5 (more precisely we obtain 
40 % at x = 0.6 and 95 % at x = 0.75). 



Figure 7. The exact 2-loop coefficient function C% s and reconstructions (approxs.) based on 
the lowest 5 moments N=2,4,6,8,10 and the sets of distributions of Eq. (|52|). The left plot shows the 
reconstruction for intermediate x. The right plot shows the quality of the reconstructed coefficient 
function as a distribution in the whole x-region, i.e. when it is convoluted with a realistic parton 
distribution A = u + u — d — d. The quark distributions correspond to the MRS (A) set [38] at its 
parametrization scale Qq = 4 GeV 2 , rif = 4. 

The present status of the effects of the higher order contributions to CV; s is illustrated in figure 8. 
The 3-loop contribution is obtained by matching the following 4 linear combinations of distributions 

£,ns,(3) 



l,x 2 ,5(l - x), 
log(x),l,x 2 ,5(l -x), 
log(x),l,x 2 ,5(l -x), 
log 2 (x),l,x,5(l -x), 
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,,53) 

to the 5 available moments (N=2,4,6,8,10), where the numerical coefficients of the terms [log (1 — 
x)/{\ — x)} + and [log 4 (l — x)/(l — x)] + are taken from Eq. (|5T|). Please notice that in the variation 
of simple distributions in the 4 linear combinations above we changed even the type of the expected 
smail-x singularities to obtain a conservative estimate of the accuracy of the reconstructions. 
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Figure 8. Higher order corrections to Fi due to the non-singlet coefficient function 
Cj? s ■ The different curves show the effects on F2 of higher order corrections to 
via the relation \F^{x, Q 2 ) - ~F$ n (x, Q 2 ) = C% s ® \{u + u-d-d), when both the 
quark distributions and a s are kept constant. The spread of the different NNNLO 
curves indicates the uncertainty of the 3-loop contribution (due to variations of the 
distributions as it is given in Eq. (|53|)). We used the parton distributions in the MS 
scheme from the MRS (A) global fit |38] (in the NL order) at its parametrization scale 
Q 2 = 4 GeV 2 and a s = .261, nj = 4. (LO = leading order, NLO = next-to-leading 
order, etc.) 
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For the non-singlet combination of the nucleon structure functions F2 P (x, Q 2 ) — -F| n (x, Q 2 ) that 
is presented in figure 8, C 2 s is the only relevant coefficient function. The different curves show the 
effects on F2 of higher order corrections to C 2 S when both the quark distributions and a s are kept 
constant. Although the change in F2 is indicated when the order of approximation for C2 changes 
(LO indicates C 2 ns ' (0) , NLO indicates C 2 S ' {0) + a s C 2 s ' {1) , etc.), it is in fact the quark distributions 
that may have to be modified to keep F2 in agreement with the experiments over a large (x,Q 2 ) 
range when higher order corrections are globally included. We noticed that also at 3-loops the two 
leading soft gluon contributions are to a large extent suppressed by 'subleading' contributions in 
the coefficient functions. 

At a relatively low energy scale of Q 2 = 4 GeV 2 of figure 8 we find that the 3-loop contribu- 
tion a^C 2 s '^ gives a sizable correction to the coefficient function C 2 S (about 1/3 of the 2-loop 

contribution a 2 C 2 s '^ in the x region 0.2 < x < 0.9). Furthermore, from figure 8 one can observe 
apparent convergence of the QCD perturbation theory up to and including the NNNL order. 

6 The analytic results of the 3-loop calculation 

Before we present the results for the anomalous dimensions and coefficient functions we should ex- 
plain our conventions for the different flavour factors that appear in the present 3-loop calculation. 



FC 2 FC U FC 02 FC$ FC 9 n 

Figure 9. Examples for diagrams in the flavour classes FC2, FCn, FC02, FC 2 , FC U . 

Diagrams in the flavour class FC2 (where both photons are attached to the quark line of the 
external quark legs) have a SU(n/) flavour factor Q 2 where the matrix Qf is the quark charge 

matrix from the electromagnetic current (Qf =diag(|,— |,— |,- • •) and tr(Qf) = Z)/=i e f where ej 
is the electromagnetic charge of a quark with the corresponding flavour). Diagrams in the flavour 
class FC\\ (where one photon is attached to a closed quark loop and the other one to the quark 
line of the external quark legs) have a flavour matrix tr^Qj^j Qf and diagrams in the flavour class 
FCq2 (flavour singlet diagrams, where both photon legs are attached to the same internal quark 
loop) have a flavour matrix tr^Q/) 1- The diagrams with external gluons are split up into two 
flavour classes: the class FC 2 has a flavour factor tr(Q 2 ) and the class FC n has a flavour factor 
(tr^)) 2 . 

To project out the non-singlet contributions one should take the flavour trace of the diagrams 
with the generators X a of SU(nj) (due to the diagonal form of Qf only the diagonal generators 
are relevant). This means that in the non-singlet projection only the diagrams from the flavour 
classes FC2 and FCn contribute. In this way the diagrams of the class FC2 receive a flavour factor 
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tv(Qf X a ), and the diagrams of the class FC\\ receive a flavour factor tr(Q /)tr(Q f\ a ) ■ Since the 
ratio of these flavour factors does not depend on the number a of a diagonal generator X a 



tr(Q/A») 
tr(Q/ 2 A Q ) 

it is possible to factorize the a-independent coefficient functions C£ S N as follows 



(54) 
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(55) 



where the non-singlet C™ N is independent of a and the standard normalization is C% S N = l + 0(a s ). 

To project out the flavour singlet contributions for diagrams with external quark legs one should 
take the trace with the unit flavour matrix. 

The explicit values for the flavour factors /Z2, fhi, fh2, fl^ fW (f° r the flavour classes FC2, 
FCn, FC02, FC2 and FC\ X respectively) that give the standard normalization of the coefficient 
functions for the non-singlet and singlet cases are given in table 2. Please notice that singlet 
flavour factors are chosen such that they reduce to unity if one substitutes for Q f the unit matrix 
Qf = diag(l, 1, 1, • • •). Please notice that the factor 3 in the non-singlet flavour factor flu originates 
from Eq. 



fh 


fhl 


fk2 
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fi 9 n 


non-singlet f 











singlet f 




1 


1 









Table 2. Values for the flavour factors ffo, fhi, flo2, /^2> f^ii- 



Below we present the results for the 3-loop anomalous dimensions and the 3-loop coefficient 
functions. Please notice that the results contain both the singlet and non-singlet contributions. To 
obtain the specific non-singlet or singlet case one should substitute the flavour factors from table 
2. In the results we have already put fli = fl\ = 1- The reader may notice an overall factor 2 
difference between our results for anomalous dimensions and the results in Ref. (|| which originates 
from our definition of the renormalization scale derivative d/(dln(fi 2 )) (see Eqs. fli~4| ) and (|i~5|)) 
versus the derivative d/(dln(fj 1 )) in Ref. [||. In the results for the coefficient functions we have 
put /j, 2 = Q 2 (in accordance with the evolution equation Eq.(I7)). The Riemann zeta function is 
written as Q n . 
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+n f ^F ( ~ 1687817250 ) + n f C A{~ 



1890 



a s (21.2666666667 + 0.6666666667^) 
+a 2 s (219.6240987654 - 24.699264321611/) 

+a 3 s (3609.3541896322 - 673.9430658122n/ - 11.2013383657^) 
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+a 2 s fl 2n f C F (-^§) 
+a 2 s n f C F (-4) 

_l„2^» r> ^ 3677 128 a \ 

, „2^2 I 4189 i 96 a \ 
+a s C F (- w + — (, 3 J 

+af /Z nn/ C F C A (-if + i^Cs + 32C 5 ) 

_L„3f; „ ^2 /496 1216 a 256 a \ 

_L„3f; „ r r ■ f H77679 1376 a 32 a 64a \ 
+a s /t02"-/^FL/ J 4(, 21870 45~^3 ~ — U ~ —(,5; 

_i_a,3 f 7 „ r»2 r 28249 17296 a ,32 a ,128an 
+a s /<02n/^F( _ ^430" ~ ToS - ^ + T^4 + — (,5 J 
_i_«3f/ „2^ f 77623 i 3232 a \ 

j_„3„ r« r< r 998153 , 17432 a 32 a , 80a \ 
+a s n / G F G j4 (-T0935- + ToB - ^ - T C4 + T^J 
i„3„ r»2 / 341701 1352 a ,32a\ 

1 „3^2/^_^ 7814 , 64 a \ 

_i_^3r» r<2 ( 3667498 46160 a , 32 a , oqra \ 
+0^0^0^(^0335 81-^3 + -3-C4 + 296C5J 

i.3/t2^ / 1013578 1 30776 a ooa 1552 a \ 
+a s O F C,A(. J0935- + -45- C,3 - 32C.4 3— C5j 

j_„3r<3( 201577 50864 a ,64 a ,416an 

+a s C F { — 7290 ioB - ^ + — U + —(,5) 

l + a s (0.4444444444) 

+a2(-2.1893004115n / // 02 + 17.6937658911 - 5.3333333333^) 
+a 3 s (-24.0920133486n//Zii - 79.044861424/Z 2n/ + 3.3255044776/Z ( 
+442.7409692714 - 165.1971095394n/ + 6.0302724150n|) 



+a s C F (^) 

+a 2 s fl 02 n f C F ( ?q%wq ) 

1 „1„ r< l 1376021 \ 

+a s n f C F {- 10800Q J 

1 „2ri n . ( 4649639 269 a \ 
+a s C F C A (- ssm 5-C3J 

-i- n 2 r< 2 ( 142624259 , 224 a \ 
-r«s u M 6480000 5 W 

+a'VZnn / C F C A (-lf + f C 3 + I6C5) 

_L„3f; „ ^2 /1232 608 a 128a \ 
+ a s/ < ll n / Cy F(.T5 15"^ 3 — S.5/ 

_La,3/7 „ r< r 1 1 20953881467 10438 a 121a \ 

+aJlo 2 n f C F CA{ i45 8 000000 ~ TT25"^3 ~ -75 U) 

1 „3fi „ a-»2 / 107812550749 161386 a ,121m 

-TU s j LQ2ii.f<^ F { 1530900OOOO 70875 ^3 T 75 

,„3f7 „2^ /• 16205516 , 1204 a \ 
+a s // 2n / O F (- 6834375 + 2025<,3j 

_l„3„ r< ^ f 36939762439 , 365878 a 314 a \ 
+a s n / F Oyl( 122472000 + "2835" ^3 ~ 
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i „ 3 « r 2 ( 283057919183 _ 222433 a , 314 a ^ 
' s f F v 3061800000 4725 ^ 15 W 

, „3„ 2 r <( 136926751 , 628 a \ 
■+a s n f U F { 87 4 80 oo + 405 ^ J 

_L^,3r< nl 1 26155194283 16471321 a , 1439 r _t_ A7nf \ 
+a s C F C A { 24494400 UTTb - ^ + W^ 4 + 4 ' U( -5j 

i„3/-(2 /-( / 1209798320299 , 21020311 a 1439 a 124 a \ 
-1-u s ^ F ^A\ 12247200000 47250 ^3 50 ^ 4 3 W 

I n 3^<3 / 4987373708399 , 119936813 a , 1439 a 1624 a \ 
-\-a s ly F [ 40824000000 283500 ^ 3 t 75 S4 3 S5J 

l + a s (6.0666666667) 

+a2(o.4858308642n / // 2 + 142.3434719201 - 16.9879135802^) 
+a^-18.2188461805n / /Zn + 16.6483484853n//Z 2 - 2.2086306890n|/Z 2 
+4169.2678883092 - 901.2351625706n / + 23.3550392440n|) 



j-nlfl r, r I 47344331 \ 

+a s Jh2n f C F ( 129654000 J 

1 „2„ n ( 1167429869\ 
+a s nfl^ F { 55566OOO > 
_L„2^< r> 1 2962167319 348a \ 

+a s C F C A { 22226400 5~W 

J-^2^2 / 44524614749 , 2036 a \ 
~r u s u M 3889620000 ~ l ~ 35 W 

1 „3f| „ r r< r A 313979 , 1028 a _i_ 320a\ 
+a s /ill^/^F^Al( 8400" + W^3 + — C5j 

,„3// „ ^2 /313979 8224 a 2560 a \ 
+a s /^iin / C F (^ T 5 - - W C3 - — C5j 

_i_^3fi „ r> C< 1 42613645803313 907166 a 484 a \ 
+a s /f 2n/UirOA(. 5145967260000 ~~ 231525 ^3 ~ 735^ 
-l-« 3 f/ „ r< 2 ( 44744425895431 1313258 a , 484 a \ 
-T^gj l02"-f^ F { 5145967260000 694575 ^3 ~r 735U; 
,„3f7 „2^ / 18418888183 , 3448 a \ 

+a s // 2n / G F (- 13613 67oooo + ISSZsW 

r* 109033584779183 , 358312 a _ 2836 a \ 
-TU s llfi^ F i^A\ 210039480000 ~r 2025 ^ 3 105 W 
i_„3„ ^2 / 11225895261923 417224 a , 2836 a \ 
-TU s ilf^ F { 51051262500 6615 ^ 3 < 10 5 S4J 

, n 3 T1 2 r « / 11510197169 , 5672 a \ 
TU s ii^ F y 388962000 2835^3; 

1 3^ r<2 ( 492030789922433 163363328 a , 34931 a , 4436 a \ 
-TU s K-jF^ A\ 280052640000 99225 ^ 3 ' 3675 S4 t 7 «i 

i„3^2^ /' 86882592642607 , 174615691 a 34931 a 1116 a \ 
~T~ u s^F^M 204205050000 385875 ^ 3 1225 t > 4 7 

I „3/^3 /" 4962846338797159 , 2402912488 a , 69862 a 4168 a \ 
-ru s ^ F { 257298363000000 ^ 3472875 ^ 3 "T 367 5 <i4 7 W 

l + a s (ll. 1767195767) 

+a;;(o.4868787285n / /Zo2 + 302.3987349886 - 28.0130504025^) 

+af (-16.1427176094n / /Zn + 24.1177881298n / // 02 - 1.5254891426nJ/Z 2 

+10069.6308450937 - 1816.3229292502^ + 42.6627311588n^) 



35 



_i_«2/7 „ ^» ^ 83382717493 
+a s fl 2n f C F ( 288054144000 J 

4-« 2 ti r< f 682775147983 \ 
-TU s lbf^F\ 24004512000 > 

, 2^ /-< / 3348220230689 _ 17071 a \ 
-\-U s i^F^A{ 19203609600 210 W 

, n 2^2 ( 203744735078509 , 2396 a \ 
-t-u s ^ F { 20163790080000 ~r 35 W 
_i_^3f/ „ r< C 1 ( 476784977 , 36443 /■ , 200 a \ 
+aJl U n f C F C A ( 8 93o 2 5o + -9450"^3 + — C5j 

,„3 f / „ n2 1 1907139908 145772 a 1600 ^ \ 
+aJl u n f C F { 13395375 - T4T75-C3 - — CsJ 

^n 3 f/ r, r r ( 91829623403140453 23440003 a 1369 a \ 
-t-u s j t02"-/^F^/ll 16 332669964800000 10716300^3 3780 t > 4 / 
,„3f7 „ ^2 /• 10181295454316382631 174416989 a , 1369 a n 
rll sJ <-02 ll -f^F\ 1257615587289600000 117879300^3 "T" 3780^4/ 
_i_„3f/ „2^ f 16959835520839 , 2017 a- \ 
-1-U-sJ l>02'l'f<^F{ 19443654720000 25515^3; 

4_„3„ ^ r ( 4344582726190968061 , 562955509 a _ 9883 a \ 
^u s a/o F^A{ 5988645653760000 2619540 ^3 315 <,4J 

,„3„ ri2 1 119458136235341463641 41139607 a , 9883 a \ 
-TU s llf^ F { 314403896822400000 523908 ^ 3 < 315 S4J 

1 n 3 r ,2 < o ( 1187935000932371 , 19766 a \ 
-TUgllf^FK. 27221116608000 8505 > 3 / 

1 3/- (A,2 / 327176399149038763 _ 1614116849 a , 2510407 a , 47044 a \ 
-ra s L/i?(_/ j4 ^ 136105583040000 793800 ^ 3 ~ r 264600 ^ 4 63 W 

4_„3/o2/-f / 5717961692522254249661 , 339247197811 a 2510407 a 1982 a \ 
-TU s ^ f ^a\ 5030462349158400000 " l " 916839000 ^3 8 8200 ^ 9 W 

I -3ri3/ 14319995334117562552891 , 33131600990227 a , 2510407 a _ aaAf \ 
^ u s^F\ 352132364441088000000 ~r" 33006204000 ^ 3 < 1323OO ^ 4 uu ^S5j 

l + a s (15.5298941799) 

W s (0.3859585393/1/ fl 02 + 470.8074190065 - 37.9248227990n/) 
+a;?(-15.0920382729n//Zii + 22.3320193843n/// 02 - 1.0363081217n}/Z 02 
+ 17162.3724471532 - 2787.2976921073n/ + 61.9117799688raj) 

a s n/(-i) 
W s n f C F (-^§ + fCs) 
W s n f C A (^ -2C 3 ) 
+a3//f 1 n^(-4 + ^C3-fC 5 ) 
+af/Zf 1 nfC A (|-fC3 + 8C5) 
+a3n / C7^(-^ + ^C3 - f C4 + 40C 5 ) 

i„3„ /-<2 / 28403 1 4148 a _i_16a 
+a s n/G F ( ilI g - + W C3 + -3-C4 

i„3„ ,^2/ 3444493 , 1828 a , 1 9 a , 10M 

+a s n/C A ( — Efm - + W C3 + 12C4 + 12C5J 

, „3„2 / a, 23291 104 a \ 
1 „3„ 2^1/130219 1 622 a \ 

+a s n f C A {^i S jo + 405C3J 
a s (— 0.5n/) 

+a2(-8.9183389614n/) 

+o2(-0.9007972776nj/if 1 - 130.7340963277n/ + 29.3793351486™}) 
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c 2 G 4 = w(-ii) 

i „2„ /-« / 6410867 i 18 a \ 

+a s n f C F (- 720000 + -C3J 

, _2„ /o / 1812647 9a \ 

+a s n/GA(- 3240000 - 5C3J 
+a2/Zf 1 nJC F (-^-^Cs + ^C6) 
W s fl 9 n n}C A (^ + if^Cs - H6C5) 

_i_«3„ r< r 1571010034879 25747 a 969 a , ooa \ 
+a s n / F L j4 ( 61236000000 630~« ~ "g^U + 8A5J 

i_„3„ ^<2 / 4566936087251 , 18867971 a , 261 a 292 a n 
-ra s ibf<^ F y 122472000000 " r 283500 S3 " r 50 ^4 3 S5J 
4_„3„ ^<2 / 77953723693 , 1312997 a 1 354 a _ 28 a \ 
~ ru s n f'^A\ 729000000 " r 20250 ^3 t 35 M 3 S5; 

1 „3„V / 790842890021 _ 595883 a ^ 
" ru s"'/ F ^ 61236000000 141750 "»3 J 
„3^2 n .1 109146757 , _59_/- \ 

+a s n/UA^ 9112500 + 2025S3J 
= a s (-0.7388888889n/) 

+a^-14.2715869197n/) 

+0^-1.6118165124^/1?! - 346.461275618471/ + 46.5201756441n|) 

C^6 = a s n /( _ 2lSo~) 

1 „2„ r» / 12660217 1 20 a \ 

+a s n f C F {- 1296540 + — C3J 

1 „2„ /o / 16794343 10a \ 
+a s n/OAi- 7779240 ~ T 

i„3f79„2^ / 1330183 1673 a i 16400 a \ 

+aJl u n f C F { mfl — C3 + -2T" C5j 

, „3 tl9 ^2 n.l 1330183 , 1673 a 2050 a \ 
+a 8 /'ll n /^l 24192 + "g - ^3 — C5j 

_l«3^ r» / 395582364061093 4535294 a 11833 a i 1040 a \ 
+a s n f L, F L,A{- 13722579360000 ~~ 231525 S3 735 _( ' 4 + ~2~rs5j 

1 3„ ri2 ( 986797608696253 1 65099719 a , 620 a 440 a \ 
-ru s iif<^ F { 13722579360000 ~r 1389150 S3 " r 147 S4 7 S5J 

i„3„ a-<2 / 230534380633807 , 43271551 a , 2911 a _ 80 a \ 
T«s' l / u J 4l 1646709523200 694575 S3 "r 2 45 S4 2lS5j 

1 „3 r) 2-^ / 8883735865961 _ 485917 a \ 
-r«s"'/ Ly - F V 571774140000 138915 S3 J 

1 „3„2/-( / 4519688643989 , 143 a \ 
~T~ a s n f^A\ 294055272000 7938S3J 

= a s (-0.7051587302n/) 
+a2(-20.0684982842n/) 

+af(-1.4960369382raj/i? 1 - 715.0372438398n/ + 61.2854509604nj) 
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n n ( 16231 'I 
u s n f\ 25200' 

i „2„ r» f 23994462871 , 7 a \ 
+a s n/G F (- 2351462400 + 3«J 

i 2 /-» I 89216988167 _ 7/ \ 
TUgllf^AK 29393280000 6^ 3 ' 

, r ,3f]9 „2r< ( 692580563 2357104 / , 15232 r \ 
+a s Jl u n f U F { 1913625 2025 - ^3 H — <,5j 

i „3 -f;3 „2/-» . / 692580563 , 294638 / 1904 / a 

+a J /f 11 n / GA( 510300Q + -675-C3 — Cs) 

,„3„ /-« /o / 13175786598173293133 _ 2254405091 / _ 515201 / , 289 / \ 
-TU s iif^F^A\ 359318739225600000 336798000 ^3 37800 ^ 4 ~r 9 W 
,„3„ z-,2 / 954577609184015777449 , 154197212633 / , 749 / 130/ \ 
-TU s ibf<^ F \ 10060924698316800000 " r 4715172000 « + 216 1 ' 4 3 W 
,„3„ z-,2 / 3720425499755558929 , 12414318401 / , 64021 / _ 10/ \ 
"""^"'/'-'Al 23332385664000000 " r 214326000 ^ 3 " r 6300 ^ 4 9 
4_„3„2/-» / 4082506317130294583 _ 187698097 / \ 
-r"s"'/ Lyir V 239545826150400000 67359600 ^3 J 
,„3 T ,2/-, / 3281715651209153 _ 4591 / 
-ru, s /tj<_yAV 194436547200000 72900 t >3 

0.6440873016n/ 



■(■ 



+o£(-23.1787352382n/ 

+0^-1.2864009150^//^ - 996.503870949671/ + 68.6646730444n|) 

a s C F (l) 
+a 2 s fl 02 n f C F (-§) 
W s n f C F {-§) 

W s c F c A (^ - f c 3 ) 

,„2/-y2 ( 1906 1 64/ \ 
+CL S C F { jgg- + — C3J 

+a3Am / C F C A (-^ + if^Cs - 192C 5 ) 
+«fAm / C2 (^-^3 + 512(5) 
+af/fen / C F C A (-^f + 3jj*Cs - Cs) 

W s flo2n}C F (^ + ff go 

_i_„3„ /o / t / 204548 1568 / , 320/ \ 

+a,n / G F Gi( — ^ 45-C3 + — Cs) 

i„3„ ^2/25534 2848 / \ 

+a a n/C F (- a5 g 45-C3J 

1 „3„2^ /2168\ 

, „3^ ^2/548668 3680 / , 99/1 / \ 
+a s G J FG A (^ 2 - T 5 — (,3 + 224C, 5 j 

_i_„3/-y2 /~t 1 41536 1 73504 / 5248 / \ 
+a s C F CA{ 405- + V3 3— <,5j 

i„3r»3 / 232798 39424 / , 3584 / \ 

+a s O F ( — 1215 45-^3 + — Cs) 

'1.7777777778") 



+a2f-3.9506172840n///o2 + 56.7553015166 - 4.5432098765n/ 
+af(-7.7366982885n / //n - 213.9253075658n/// 2 + 17.9132652795n 2 ,// 02 
+2544.5980873744 - 421.6908884762n/ + 11.8957475995™^) 
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a s C F ( 4 



586 



+ajfl 02 n f C F (-^) 
W s n f C F (-%) 

^r,2r< n / 4763 48a \ 
i r ,2r<2 ( 19967 i 96 a \ 

+a3/Z 11 n / C F C A (ifP + ^Cs - 192C 5 ) 

_i_«3f/ „ r 12 ( 49436 6848 a , rio/-\ 
_i_«3f7 r> r> l 5074537 824 a \ 

+aJlo2nfC F L A {- 405Q00 - 675C3J 

_i_^3f/ „ 7715012 , 29488 a \ 

+a s /<02n/O F (- 1063125 + ^725-C3j 

1 „3f/ „2^ /■ 106981 -> 
_l_«3„ n TJ 14259893 , 55904 a \ 

+a s n f C F C A { 94500- + -945-C3) 

,„3„ r*2 ( 258828431 _ 33344 a \ 
-TU s itf<^ F { 2835000 315 W 

+afn^ F (ifi) 

j-^Zr 1 n2 r 171354151 1889656 a , 00 a \ 
+a s G jF A ( 283 5oo 4725- C3 + 32(^5 J 

, .3/^2/1 f 263016326 , 1070306 r 1 09/- \ 
+a s C F G A ( 354375 + 1575 C3 + 32C 5 ) 

1 „3^3 / 582157141 , 57356 ^ 1QV \ 

a s (l. 0666666667) 

+a2f-0.6945185185n//Z 2 + 47.9939893120 - 3.4133333333n/) 
+a3(-5.0588695123n///n - 55.553045597171/7/02 + 2.3480054870n}/Zo2 
+2523.7390200791 - 383.0520013416ra/ + 10.8889547325™}) 



a s C F ( 4 



7> 

14614 



+a s fl Q2 rif C F (- jytiE ) 

, _2„ r< / 1486 \ 

+a s n/GF(-735-J 

_i_«2 ( ay a~y / 172106 48 a \ 

+a s C jP GA( TJQ25" ~ T«J 

_i_„2ri2 ( 257318 1 96 a \ 
+a s C F (^- 25725 + — (, 3 J 

_L„3f; „ r 1 n . ( 5869993 1 944 a 960 a \ 

+aJlun f C F CA{ 6 3ooo + ^rC3 - —Qb) 

_i_^,3f7 „ ^2/ 5869993 7552 a , 2560 a \ 
+aJl n n f C F { 23025 gg-Cs + — C5j 

_i_«3f/ „ r< r,( 41368883483 7384 a \ 
+a s Sl 2nfC F C A {- 8 i 68 202000 - TT025«J 

_l„3 f 7 „ ^2 / 193467963407 , 416 /■ \ 
+a s Jh2n f C F {- 6 i 26 i5i5ooo + 175^3) 

, „3 „2/-f ( 27905687 \ 
+a s Jl 02 n f C F { 48620250-' 

_i_^.3^, r< r< f 927531079 , 88916 a \ 
+a s n/G F GA( 694 5750 + "lW^J 

_i_„3„ ,ay2 /■ 66489992539 _ 31312 a ^ 
TU s iif^ F ^ 875164500 315 W 

I n 3 T7 2 rrr( - 986872 x 
+ a s n /°^l 138915 
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j-„3r< nl 3186722339 4977821 / , 1040 a \ 

+a s u F u A { 6174000 rTo25~^ 3 + 

j_ n 3r2r,( 4333877446411 , 2071508 a A o n /- \ 
+a s C F C A ( 8751645000 + 2205 « ~ 480<, 5 J 

-4-n 3 r< 3 ( 90265366481 _ 536548 a , 2560 a \ 
~r u s^F\ 2268945000 2205 ^ 3 7 t > 5 ^ 

a s (o.7619047619) 

+ a 2(-0.2524824533n / // 2 + 40.9961975991 - 2.6956916100n/ 
+a 3 s (-3.7056125257^/^1 - 24.0132253893n//Z 02 + 0.7652692585n^// 02 
+2368.1937754336 - 340.0691069253n/ + 9.4721904282n 2 / 



a s C F ( 4 



+a 2 Jl 02 n f C F (-^) 

_i_«2.ay / 14741729 16/- \ 
+a s C jF CA( H90700 ~ T^j 
j_„2r<2 ( 21694349 , 32/- \ 
+a s°F^ 3572100 ~*~ T«J 

_l^3 f 7 „ r< r> 1 35555777437 , 170756 /- 320 a \ 

+a s //iin / G F GA(^ooo940or + T725-<>3 - — C5J 

_i_^,3f7 r<2 f 35555777437 1366048 a , 2560 a \ 
+a s //nn/G F ( 1 87535250 14175 - ^ 3 + — C5j 

_l„3^7 „ n r< ( 4264597461073 12892 a \ 
+aJlo2nfC F CA{- 1620304560000 ~ 33075^) 
4-n3f7 „ r<2 I 55595717822203 , 1327432 a \ 
-r«sJ 31190862780000 ^ 1091475 *»3 J 
, „3 „2^y / 3487324703 \ 
-ta s Jl 02 n f U F { 13502538OOO > 

ri ri ( 23153083641529 , 18459136 /-_\ 
f^F^A{ 198037224000 " + " 



S »/^^A\ 198037224000 1 363825 
3^ r<2 ( 11844644404289 914992 a \ 



198037224000 10395 W 
I „3 T? 2 ( o / 1435876 \ 
+ a s n /°^V 229635 ^ 

i flp fi2 ( 7653142193467 _ 9508318 a , 2240 a \ 
-TU s i^ F ^ A { i 8 oo3384000 19845 ^ 3 T 9 

_I_„3/a,2/-» / 7667007621800089 , 2476882549 a 2720 a 
-ru s ^ F ^A{ 27725211360000 " r 2182950 ^ 3 S5 
_i_„3/^3 / 86167166469418457 _ 111668693 a , 7360 a \ 
-TU s i^ F { 499053804480000 218295 ^3 ' 9 ^5) 

a s (0.5925925926) 

+^(-0.1217397796^/^02 + 35.8766440564 - 2.2314716833^' 
+a;?(-2.9137025628n / /Zii - 12.9718526710n / /Z 02 + 0.3443623910n^/Z 02 
+2215.2108750618 - 305.4730328944n/ + 8.3371495344n} 
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Cg 2 = a 8 n f {\) 

+a*n,C F {-ffi - f Cs) 

+a 2 n/C . A( m ) 

Wjl 9 n n}C A (-l - fks + H2C 5 ) 

^3„ r» r< I 71657 248/- i 80/- \ 

i„3„ ^2^ 51283 i 928^ 160 a \ 
+a s n f L F { Tm - + ^-C3 - — C5j 

i „3„ /o2 ( 235283 148/ , 64/ n 

+a^njC F (m + WC3) 

W + &Cs) 
= a s (0.6666666667n/) 

+a^l2.9477670897n/) 

+0^-0.3889396640^/^^ + 407.280631959671/ - 20.2395974792n^) 

= «<«/(£) 

+a^/C F (-ii) 

1 „2„ / 19229 \ 

+a3 / /f 1 nfCH-^ + C 3 ) 
+af/Zf 1 n 2 ^(^ - *f Cs) 

_i_„3„ r> r> { 261555163 272858 / , 100/ \ 

+a s n f C F L A ( — isgoooo 4725~C3 + 1^5 J 

_L„3„ nl 1 912595079 1 225908 /- 190/ \ 

+a s n/C F ( 11340000 + T725-C3 - 128C5J 

_i_„3„ r 1 " 2 ( 172349813 , 5036 / 09/ \ 
+a s n f L A { 12150Q0 + W C3 - 3^C 5 J 

1 „3^2/^_/ 4921307 32 /- \ 
+a s n/OFl, 680400 - 105S3 J 

1 (1 3„2 / o / 445666 8 / \ 

= a s (0.2666666667n/) 
+a2(l3.8165925926n/) 

+af (-0.3984298404n}/Zf 1 + 767.7125420910n/ - 36.7841923177n^) 

Cl,6 = W7) 

1 „2„ /o / 137761 \ 

j_„2„ r 1 < 503977 \ 

_i_„3 f/9„ 2^< / 7994213 , 61778 / 1760/ \ 

-\-a s ji u rifU F { 340200 + "STB - ^ — <,5j 

x«3f/9 „2p / 7994213 30889/ _j_ 660/\ 
+a s jL 11 n^A\- 907200 420~^3 + — t,5j 

_l^3^, n n ( 1620883880209 51181/ , 480 a \ 

+a s n f C F C A { 1 6 33 6 404000 - ^20r^3 + ~ 

_l_„3„ r<2 ( 101548438631 , 43574 / 480/ \ 
T«s' t / Ly Fl 2042050500 2205 ^ 3 7 <,5J 
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i 6 ri'l f 38038WU2Y6Y , IU»U/ _ 
-ra s Hf^ A y 3267280800 ^ 4410 ^ 7 ^ 5 

, n 3„2^ / 60190474091 _ 32 a \ 
TU s"/ F v 12252303000 63 W 

, „3„2/-< / 847812559 2 a \ 
+a s fl f L, A {- 77792400 ~~ 7^) 

= a s (0.1428571429n/) 
+^(10.260953640771/ 
+af (-0.3055276256n^ /Z?! + 694.5092121136ti/ - 27.9895081033ti}) 

C L,8 = a snf(^) 

i „2„ r< l 51097 a 

i „2„ ^ ,/ 7712869 \ 
+a s n/GA( 2551500^ 

< r ,3fl9„2r< ^ 3665714041 , 617672 a 1664a \ 
+a s Jl n n f L, F { 35721000 + 1575 <>3 —^5) 
_l_«3 f/3 „2r< , ( 3665714041 77209 a , ono/ \ 

+a,// 11 n / C j4 (- 952560 oo 525-« + 208C 5 ) 

4-n 3 « r r ( 520855237960033 _ 6119609 a , 128 a \ 
tUj/ijL-FUi^ 7129340064000 519750 W 1 " 3 ^> 

i 3 r~i2 ( 2384408424295187 , 7723411 a _ 128 a n 
" l_u, s"'/ Ly -Fl 71293400640000 779625 ^ 3 ^> 

I „3„ ,ay2 / 27404278602137 , 20438 a _ 32 a \ 
" l_u 's' t / < - y Al 289340100000 ~r 23625^3 3 ** 5 ) 
_l„3„2 ( ay / 124374980290567 _ 608 a a 
-TU s ilfi^F{ 35646700320000 1485 S3 J 

i „3^2^ . /- 11324757281 8 a\ 
"1" a s ri 7 U M 1377810000 45^3; 

= a s (o.0888888889n/) 
+^(7.733545104271/) 

+a3(-0.2322211886n^/Z? 1 + 592.330797209871/ - 2 1.3033368 11 In)) 



In addition to the results given above for the moments N=2,4,6,8, we present the following results 
for the non-singlet contributions to the 10th moment. 



7fo s 



a s C F ( 



12055 - 
1386 - 



r< n 1 19524247733 \ , r 1 ^ t 9579051036701 
^F^A{ 523908000 ) + U F\ 



(J p (j1 (94091568579766453 



435681892800000 

I ~r<3 ( 2207711300808736405687 i 151796299 a a , p /-i t 

"r^M 127866318149354400000 ^ 4002075 s3 ; < uf<^F<^A\ 389001690000 
+r> C 2 ( 75522073210471127 , 48220 a \ , 2 fi t 27995901056887 
/ 1230075210672000 ~r 693 ^> ^ n f K - y F\ 11981252052000 



_l_ n r ( 2451995507 \ 
FV 1331250228000^ ~r n f^F\ 288149400 / 

151796299 a_ A i r<2 n . t 16389982059548833 



8004150 £3) + C f Ca{ 465937579800000 
_|_ 151796299 a„ A , ^,n„n.( 900777312740o 



151796299 

2668050 
_ 48220 a 
693 s: 



GO 



a s 11.5969215969 



W s (136.2741775192 - 11.3459453418ti/) 

W s (2379.9199516669 - 387.6422968015ti/ - 3.1155231451n}) 
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Cns -i 
2,10 — 1 



2006299 - 



138600 

/ 56: 

2n n ( 6124093193824187 104674, 



+a s C F ( 

i 2 fi i 561457267429757 \ 
-ru s iif<^F\ 15975002736000 / 
\ n 2(-i (~i (6 

-ru s ^F^A{ 29045459520000 1155 ScW 
i 2/i2 ( 558708799987324013 , 88798 a \ 
+u s L 'i f i 14760902528064000 1155 S3J 
,„3fi „ r> ( 3753913187503 162776 /- , 896 a \ 
+a s /tnn/0 F OAi 58677696000 202125 S3 + TT S5j 

4-n 3 f7 n r 2 ( 3753913187503 , 1302208 a 7168 a \ 
TU S J lullf^ F y 22004136000 " r 606375 S3 33 S5J 

4_„3„ fi fi ( 21664244926039357214987 , 10519793104 a _ 24110 a \ 
-TU s ltf^F^A\ 23550349033411200000 ~r 42567525 S3 693 S4J 
i_„3„ /o2 / 1521387460036994061010049 _ 3997754476 a , 24110 a \ 
tlt s"/°Fl 2720065313358993600000 42567525 S3 "T 693 <,4J 

_l_„3„2^ / 57084428047851551911 , 48220 a \ 
-TU s Uj^f{ 99 6 3 6 0920644320000 18711 S3; 

J-/) 3 /^ r 709221119965457939095237 14713925739913 a , 151796299 a , 190858 a \ 
-ru s ^F^A\ 235503490334112000000 6243237000 S3 T 10008300 S4 T 2 31 S5J 

i„3^2^ 1 16350009304926933389608829 , 1430215936081 a 151796299 a 22658 a \ 
+(1 s u F L/ ' 4 l 8369431733412288000000 ~r 6163195500 S3 5336100 S4 99 <,5J 

I „3f3 ( 3247779532370920623770610131 , 2182208825245282 a , 151796299 a 75212 a \ 
TU S ^ F \ 92155812816602703168000000 ~r 1622461215375 S3 T §004150 S4 99 <,5J 

1 + a s (19.3006156806) 

+a 2 s (639.2106629599 - 46.8613184159^) 

+a 3 s (- 14.4587445075^/^1! + 24953.1349702005 

-3770.1021201303n/ + 80.5209797251n|) 



,„2„ n ( 163679 \ 

^r, 2 r n . 1 89670761 48/ \ 
+a s C F CA { 8 73i 8 oo ~ TTS3J 

, „2,ay2 ( 1999510607 , 96 a > 
-ru s ^ F \ 528273900 11 S3; 

_i_^ 3 f/ „ r> C ^ 5073093424963 , 3641546 a 960 a 
+a s JtiinfU F L,A{ 88016544000 + 121275 S3 ~ — s! 
4-« 3 f/ r> r 2 ( 5073093424963 29132368 a , 2560 a 
TU s j liiitfi^py 33006204000 363825 S3 T u 1,5 

i„3„ ri (~t A ( 176183576988227323 , 55485434 a \ 
TU s nf^F^A\ 1699159381920000 1216215 S3 J 
! fi2 1 9048874326307637 _ 1174256 a \ 
i n f^F\ 190368782604000 15015 S3j 
1 n 3„2r</ 63272639 \ 



,5) 



+a 3 



's 
,3„2 



, ^ s .uf^ r v 11320155^ 

1 3/o ^2 f 2366034921481985137 95022195887 a , 3760 a \ 

i-UjUFU^^ 6796637527680000 187297110 S3 ~T n <,5j 

I n 3fi2 fi ( 323139848004267269 , 22904191 a 14240 a \ 

-TU s ^ f ^a\ 3354750574560000 17325 S3 n S5J 

1 3^3 / 887562386698645967383 _ 357031607224 a , 13440 a \ 

~ra s ^ F \ 3106213592269728000 468242775 S3 ' n S5J 

a s (0.4848484848) 

+a 2 s (32.0176594698 - 1.9085982480n/) 

+a 3 s (-2.397641694571/ fl n + 2081.2132221274 

-278.0172176870n/ + 7.4525056120n^) 
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